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We demonstrate that the Nambu-Goto string spectroscopy with massless quarks is replicated 
in highly excited states of the linear vector confinement potential. For deep radial excitations we 
observe that the Regge slope, spacing between daughter trajectories, and absolute state energies 
agree with those of the QCD string. 



I. INTRODUCTION 

We present here an interesting simple observation on 
electric potential confinement. Our conclusions concern 
the nature of potential models and their relationship to 
QCD string confinement. In the case of the time compo- 
nent vector (TCV) potential we note a remarkable tran- 
sition in the Regge spectroscopy from the leading trajec- 
tories with slopes l/4a to deep daughter slopes of 1/na. 
In fact, the exact string spectroscopy is found to be em- 
bedded in this regime. The latter slope is characteristic 
of the Nambu-Goto QCD string with one fixed end pj. 
We first discuss some general properties of the spinless 
Salpeter (SS) or square root equation. For the linear 
TCV potential we then compute the Regge structure in 
the orbitally and radially dominant regimes. 



II. SPINLESS SALPETER EQUATION 

The spinless Salpeter equation has long been attrac- 
tive as a minimal generalization of the non-relativistic 
Schrodinger equations with relativistic kinematics. Its 
main difficulty is that it involves the square root of a dif- 
ferential operator, sometimes called a pseudo-differential 
operator, which is difficult to treat analytically. Even so, 
it is easily implemented numerically. The introduction of 
a variational method allowed straightforward numer- 
ical solutions involving square roots and even transcen- 
dental functions of operators. In this section we discuss, 
though do not prove, the necessity of using the spinless 
Salpeter equation in place of a Klein-Gordon type equa- 
tion when dealing with a pure vector potential. 

Starting from the simple classical Lagrangian, 



L = -(m + S(r))V 1 - v 2 - V{r), 



(2.1) 



describing a particle of mass m moving with velocity v 
in a Lorentz scalar potential S(r) and a Lorentz time- 



component vector potential V(r), we find that the Hamil- 
tonian, 



H = ±y/p 2 + (m + S{r)) 2 + V(r), 



(2.2) 



has a sign ambiguity. This ambiguity arises from the 
elimination of the velocity in favor of the momentum 



dL 
9v 



(m + S(r)) 7 v, 



(2.3) 



where 7 = (1 — v 2 ) -1 - 72 . The positive sign in Eq. I|2.2() 
corresponds to the usual positive energy states and yields 
the correct non-relativistic limit 



H 



\p\<r. 



I_ + S (r)+V(r). 
2m 



The squared form of Eq. (|2.2|) , 

(H-V(r)) 2 =p 2 4 



(2.4) 



(2.5) 



known as the generalized Klein-Gordon (GKG) equation, 
behaves quite differently in the cases of a pure vector and 
a pure scalar potential. This difference is most easily seen 
from consideration of the s-wave turning points, where 
p 2 = 0. Eq. i|2.2[) determines the energy as a function of 
the potentials at the turning points, 



E = ±(m + S)+V. 



(2.6) 



We show this turning point condition separately for 
pure scalar and pure vector linear confining potentials in 
Figs.^andEl In Fig.^we see that for the pure scalar case 
there is no overlap of the positive and negative branches 
of Eq. I|2.6I) . In Fig. [3 on the other hand, we see that 
for any E > m in pure TCV linear confinement there are 
states on both positive and negative branches. A com- 
pletely analogous observation [3| has been made earlier 
in the context of the Dirac equation. 
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Since the GKG equation l|2.5[l is quadratic and con- 
tains both branches of Eq. I|2.2|) , it can have no normal- 
izable solutions in pure TCV confinement because there 
is tunneling to the negative energy states, which becomes 
catastrophic for small quark mass. This difficulty is the 
well-known "Klein Paradox" jj] . On the other hand, the 
GKG equation with pure scalar confinement does have 
well defined normalizable solutions because there is no 
possibility of tunneling to the negative branch. 

The spirit of the Salpeter equation || is to conserve a 
definite particle number. To this end, energy projection 
operators are employed to remove the negative energy 
spectrum and the fermion TCV Salpeter equation has 
normalizable solutions for confined quarks |3J. The spin- 
less Salpeter equation for TCV confinement achieves the 
same end more simply by explicitly including only the 
positive branch of Eq. (|2.2() . For p-wave and higher an- 
gular momentum states, we consider the Hamiltonian in 
spherical coordinates, 



H 



Pr 



J2 



(2.7) 



which we obtain from dividing the momentum into radial 
and angular parts. The wave equation to be solved is then 
the eigenvalue problem 

H\tP) = E\tP) (2.8) 

for the pseudo-differential operator given by Eq. (|2.7|) 
with the usual operator replacements. 



III. SEMI-CLASSICAL ANALYTIC 
QUANTIZATION 

We have discussed the necessity of using the SS equa- 
tion instead of the GKG equation with vector confine- 
ment. In semi-classical language, the necessity of using 
the SS equation becomes the necessity of using the turn- 
ing points given by Eq. i|2.6fl with the positive sign only. 
In this section we calculate the semi-classical spectrum of 
the SS equation with TCV confinement. We first rewrite 
the SS equation (|2.7H in terms of dimensionless variables 

ar 

X = E' 



B = — 



(3.1) 



Pr 



Pr 

E' 



In the case of a massless quark, we have 

Pr = — , (3.2) 

X 

Q = (x- xi-)(x - x 2 -)(xi + - x)(x 2 + - x), 



where 



x 1± = ^(liVl - 4^), 
x 2 ± = \ (l ± V± + 4/3) ■ 



(3.3) 



The dimensionless physical turning points resulting from 
the SS equation are x\±. The dimensionless constants 
X2± are turning points of the GKG equation, but not of 
the SS equation. 

In the regime we are investigating the angular momen- 
tum J is limited but the energy E is large. It follows 
that (3^1. To see this, we note that even for the lead- 
ing (circular orbit) Regge trajectory J/E 2 = l/4a (see 
the appendix) and hence /3i ea ding = 1/4. For fixed J and 
large E, B becomes small. In this limit the leading be- 
havior of the roots (|3.3I) is 



xi- = (3 + (3 2 + 



xi+ = 1 - & ~ P 2 ~ 



X 2 - = -P + (3 2 + . . . , x 2 + = 1 + - (3 2 + 



(3.4) 



The SS equation may be quantized semi-classically by 
the usual method 6] 



T+ j ( 1 

ar p r = 7r I n H — 



(3.5) 



In dimensionless variables l|3.1|l . the quantization condi- 
tion (|3.5|l becomes 



dxP r (x) = -g2 [ n +2 



(3.6) 



The TCV quantization is achieved by integrating 
Eq. I|3.2jl . The method is quite accurate for all states 
and becomes exact for states with many nodes in the ra- 
dial wavefunction. This is the regime we are particularly 
interested in. 

Although the integral of Eqs. (EH) 

can be ex- 
pressed in terms of elliptic integrals, it is more efficient 
to approximate the integrand first for small (3. The ap- 
proximation 



Q \J(x~ X\S)(x - X2-) 



1 



B 2 



2(1 -x) 

(3.7) 

reproduces the exact integrand to better than 0.5% 
throughout most of the region of integration, even when 
[3 = 0.1. The accuracy of this approximation is shown 
graphically in Fig. 3, where we plot both the exact inte- 
grand i,'>.2\) and the approximation (|3.7|) for (3 = 0.1. 

With the approximation l|3.7|) . the quantization inte- 
gral (|3.6|l is relatively easily evaluated. To order f3 we 
find 



(3.8) 



which immediately leads to the result 
E 2 

— = J + 2n+l. (3.9) 
ira 

Upon making the Langer correction to take into ac- 
count the centrifugal singularity 



J -> J 



(3.10) 
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where J is now the angular momentum quantum number, 
we find the final spectroscopic relation 

E 2 3 

— = J + 2n+-. (3.11) 
na 2 

The relation l|3.11|l is identical to an analytic approxi- 
mation for highly radially excited QCD string, as well as 
being an excellent approximation to its exact numerical 
solution @. 

IV. COMPARISON TO EXACT NUMERICAL 
RESULTS 

We can nail down our central result that the QCD 
string spectrum l|3.11|l is replicated in the radially dom- 
inant regime of the TCV potential by solving the 
TCV SS equation exactly numerically. The variational 
(Galerkin) method is well suited for solving eigenvalue 
equations with mixed coordinate and momentum opera- 
tors. Briefly, the method begins with a complete set of 
orthogonal states that can be Fourier transformed. The 
wavefunction is approximated as a superposition of the 
lowest N of those states. The wave equation transforms 
into an N x N matrix equation which is then diagonal- 
ized. The accuracy of the resulting eigenvalues and wave- 
functions is measured by the dependence on N and the 
dependence of the scale parameter of the basis set. De- 
tails can be found in 9]. Some further considerations of 
the accuracy of the variational method are presented in 

m- 

The TCV SS equation (|27l) with m = can thus be 
solved for a variety of angular momenta, J = 0, 1, 2, . . . 

and radial states, n = 0, 1, 2, The result is depicted 

in Fig. 4. Plotted on this figure are both the QCD string 
slope I/to (solid) and the usual leading trajectory TCV 
slope l/4a (dashed). In both cases we start the lines at 
J = 0. Although the leading trajectory (n — 0) agrees 
with the dashed prediction, the situation changes as we 
examine the higher radial excitations. At the larger n 
values (deep daughters) we observe that the solid line 
corresponding to the QCD string slope becomes accurate. 

In Table 1 we provide the exact and the string solutions 
for the s-wave radial excitations. We see now that all 
aspects of the QCD string result (|3.11|) work quite well. 
The difference between excitations is accurately two units 
of E 2 /na and the absolute values of the state energies 
correspond well to the QCD string to three significant 
figures. 

V. SUMMARY AND DISCUSSION 

On occasion, some simple results are quite unexpected. 
We did not expect to find the m — QCD string spec- 
troscopy in a potential model. To summarize our finding, 
we have examined both numerically and analytically the 
state spectroscopy for a massless quark moving in the 



time-component vector linear confinement potential. For 
states with angular momentum much less than the num- 
ber of nodes we find the spectrum is exactly that of the 
QCD string. The Regge slope, radial excitation energy, 
and the absolute values of the energy (|3.11l) are exactly 
what one expects from the QCD string. 

Simple results, even surprising ones, usually have sim- 
ple explanations. As we have recently pointed out 
the QCD string equations reduce to the spinless Salpeter 
equation with a linear time-component vector potential 
for the s-wave states. Thus there is a natural physical 
connection between the two systems. We can understand 
this connection simply. A QCD string corresponds to a 
constant chromoelectric field in the quark rest frame. In 
the limit that the quark is moving radially, the QCD 
string has no angular momentum. The total energy of 
the system becomes 

,-, r^> o arcsinuj_ v± -*a , — = t 

-C'strine = VP + m +CLT ► V Pr + m z +ar. 

(5.1) 

In a potential model the field never carries angular mo- 
mentum. The potential energy of the quark is linear in 
the distance from the origin, which is the same as the 
energy of the string as long as the string is moving ra- 
dially so that there are no relativistic corrections due to 
its transverse motion. The total energy of the quark in 
a TCV potential, 



E TC v = vV + ni 2 + ar, (5.2) 
in the limit of vanishing v± becomes 



lim String = hm i?TCV = \/vr 2 + m 2 + ar. (5.3) 

It is thus natural to expect agreement between string and 
TCV confinement in the radial dominant regime. 

The real puzzle is why the Regge slope of the QCD 
string should be the same for both circular and radial 
motions. On one hand it is well known from Nambu-Goto 
days that the Regge slope of a rotating QCD string with 
one end fixed is 1/wa. We have seen that the TCV po- 
tential Regge slope varies from 1 /4a to 1/ira as one goes 
from orbital to radial motion [see Fig. 4]. There is thus 
no obvious reason why this should not happen with the 
QCD string. Since the QCD string/TCV potential Regge 
slopes coincide at 1/ira for nearly radial motions this ex- 
plains the remarkable uniform Regge structure seen pre- 
viously from numerical solutions of the QCD string |8| ■ 
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APPENDIX A: NEARLY CIRCULAR ORBITS 

For completeness we semi-classically quantize the TCV 
equation for nearly circular orbits (J 3> n). Using the 
notation of Eqs. I|3.1H3.3[I . we consider the orbitally dom- 
inant regime. In this case the turning points x\+ and x\- 
are nearly equal. By Eq. i|3.3[) . this occurs exactly at 

1 

Xc ~ 2' 
1 

4' 



(Al) 



For nearly circular orbits we can expand in /3 and x, 
giving 



Xl± 



Q 



16 2 V 2 



(A2) 



- 2^ ~ x i-)( x i+ ~ x ) 
Using the quantization condition l|3.fi|l . we obtain 



2^2 



(A3) 



with /3 = aJ/E 2 . Solving for y/]3, squaring and dropping 
the small squared n/E 2 term, we obtain 

^ = J + V2n+-^. (A4) 

Finally, we make the Langer [7j correction J — > J + 5 
and find [g] 



£ 2 r 1 1 

— = J + V2n + - + 
4a 2 V2 



(A5) 
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TABLE I: An exact numerical solution for s-wave TCV con- 
finement compared with the analytic approximation to the 
TCV/string showing the rapid convergence of our WKB ap- 
proximation as the number of radial nodes, n, increases. 



n 


Exact (f^) 


2n+l 





1.59 


1.50 


1 


3.53 


3.50 


2 


5.52 


5.50 


3 


7.51 


7.50 


4 


9.51 


9.50 


5 


11.51 


11.50 


6 


13.51 


13.50 


7 


15.51 


15.50 


8 


17.51 


17.50 


9 


19.51 


19.50 


10 


21.51 


21.50 




FIG. 1: S-wave classical turning points for Eq. 12.61 with 
scalar confinement. 




FIG. 3: The exact Bohr-Sommerfeld integrand l:->.2l and the 
approximation 113.711 for j3 = 0.1. 




FIG. 4: Regge diagram of exact numerical solutions to 
Eq. 12.811 with linear vector confinement (dots). The solid 
lines are the QCD string trajectories with slope 1/na and the 
dashed lines have slope l/4a. We observe the transition from 
the leading TCV slope to string slope. 



